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1.0  INTRODUCTION 

The  interaction  of  plasmas  with  electromagnetic  fields  is  a  topic  of  considerable 
breadth  and  increasing  importance.  The  Air  Force  Rocket  Propulsion  Laboratory  (RPL) 
at  Edwards  Air  Force  Base  is  particularly  concerned  with  the  interaction  of  microwave 
signals  and  rocket  exhaust  plumes.  It  is  known  that  the  free  electrons  in  rocket  exhaust 
plumes  interact  with  radar  signals  and  may  result  in  observable  plume  radar  echoes.  Also, 
the  plume  can  severely  interfere  with  the  transmission  of  microwave  signals  through  the 
plume  by  attenuating  the  signal  or  by  phase  distortion.  It  is  impractical,  if  not  impossible, 
to  measure  the  radar  cross  section  (RCS)  of  all  missile  plumes  over  the  range  of  flight 
conditions  of  interest.  Therefore,  accurate  and  reliable  models  are  required  to  predict  the 
plume  RCS.  The  problem  of  predicting  plume  RCS  is  twofold.  First,  the  electrical  properties 
of  the  plume  must  be  modeled.  This  requires  complex  gas  dynamic  and  chemical  analysis 
of  interactions  of  the  external  air  mixing  with  the  exhaust  plume.  Second,  the  scattering 
of  the  microwave  signal  from  these  plumes  must  be  modeled  from  the  predicted  plume 
electrical  properties.  Both  plume  electrical  properties  models  (Ref.  1)  and  plume  radar 
scattering  models  (Ref.  2)  exist. 

A  key  parameter  in  these  models  is  the  local  value  of  electron  number  density.  It 
is  therefore  necessary  to  measure  the  local  value  of  the  electron  density  of  plumes  in 
order  to  verify  these  models.  There  are  three  basic  measuring  techniques  which  can  be 
used  to  obtain  plume  electron  density:  (1)  probes  (e.g.,  Langmuir  probes),  (2)  spectral 
emission  meaurements,  and  (3)  microwave  attenuation  or  phase  shift  measurements.  The 
severe  environment  of  rocket  exhaust  plumes  prohibits  the  use  of  probes  for  reliable 
electron  density  measurements.  For  solid  rocket  motor  plumes  the  unsteadiness  of  the 
plume,  coupled  with  solid  particulate  continuum  radiation  and  the  low  intensity  of  the 
radiation  at  the  desired  wavelength  from  the  appropriate  atoms,  makes  it  virtually 
impossible  to  use  spectral  emission  data  for  calculating  the  plume  electron  density.  This 
therefore  dictates  that  microwave  attenuation  and/or  phase-detecting  techniques  be  used 
to  measure  the  electron  density  in  solid  rocket  exhaust  plumes. 

Additional  problems  arise  when  the  electron  density  is  not  uniform  across  the  plume 
because  the  resultant  signal  represents  the  integrated  effects  of  the  interaction  of  the 
transmitted  microwave-  signal  and  the  plume  as  it  passes  through  the  plume.  It  then  becomes 
necessary  to  invert  the  measured  data  from  integrated  results  to  local  values  of  an 
absorption  eoefhcient.  This  inversion  requires  scanning  of  the  microwave  beam  across  the 
plume,  which  in  turn  generates  further  constraints  on  the  system  (e.g.,  adequate  resolution 
of  the  density  profile  requires  that  the  beam  diameter  must  be  much  smaller  than  the 
plume  diameter).  The  plume  diameter  encountered  in  most  laboratory  experiments  is 
relatively  small  (4  to  8  in.),  and  it  is  therefore  necessary  to  have  the  microwave  frequency 
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as  high  as  possible  so  that  better  focusing  techniques  can  be  applied.  These  constraints 
make  it  extremely  difficult  to  make  phase-shift  measurements  because  of  the  signal  path 
around  the  plasma  which  is  required  for  phase-shift  comparison.  At  the  higher  frequencies 
the  wave  guide  becomes  quite  sensitive  to  vibrations,  and  at  the  lower  frequencies  resolution 
is  lost.  The  method  of  obtaining  electron  density  profiles  of  rocket  exhaust  plumes 
therefore  narrows  down  to  using  microwave  attenuation  data  where  the  selection  of  the 
microwave  frequency  is  dictated  more  by  resolution  requirements-  than  by  the  plasma 
frequency. 

The  primary  purpose  of  the  project  reported  herein  was  to  examine  the  prospect 
of  obtaining  plasma  electron  number  density  profiles  from  microwave  attenuation 
measurements.  To  accomplish  this  goal  a  low  density  supersonic  argon  plasma,  produced 
by  an  arc  jet,  was  used  to  provide  axisymmetric  electron  density  profiles  from  which 
spectral  data  and  microwave  attenuation  data  could  be  obtained.  The  electron  density 
and  temperature  profiles  were  obtained  from  inversion  of  the  spectral  line  data,  using 
previously  developed  methods  (Ref.  3)  to  evaluate  the  inversion  of  the  integrated 
attenuation  data  and  the  available  theories  for  relating  the  inverted  attenuation  data  to 
electron  number  density.  The  emphasis  of  this  investigation  was  to  obtain  experimental 
microwave  attenuation  data  under  known  operating  conditions  and  to  use  existing 
relationships  to  obtain  electron  density  from  the  attenuated  data.  It  was  not  the  express 
purpose  of  this  project  to  delve  into  the  theoretical  aspects  of  waves  propagating  through 
real  plasmas,  but  rather  to  provide  good  experimental  data  for  present  or  future  evaluation 
of  wave -propagating  models.  The  results  of  the  investigation,  however,  did  point  to  the 
necessity  for  more  refined  analysis  of  the  wave  propagation. 

2.0  APPARATUS 

The  experimental  apparatus  consisted  of  an  arc  jet  which  exhausted  into  a  low  density 
test  cell,  a  spectroscopic  measurement  system,  and  a  microwave  transmitter  and  receiver 
system.  The  arrangement  of  the  apparatus  is  shown  in  Fig.  1.  Simultaneous  microwave 
attenuation  and  spectral  measurements  were  made  to  avoid  the  problems  associated  with 
trying  to  establish  repeatable  plasma  conditions.  Further,  a  square,  rather  than  cylindrical, 
quartz  test  section  was  chosen  to  prevent  lens  effects  as  the  plume  is  scanned  by  either 
the  microwave  beam  or  spectrometer.  The  microwave  system  was  mounted  on  a  scanning 
table  which  is  prthogonal  to  the  spectrometer  as  shown  in  Fig.  1.  The  experimental 
arrangement  is  further  illustrated  by  the  photograph  in  Fig.  2,  and  a  sketch  of  the  test 
cell  cross  section  is  shown  in  Fig.  3. 
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Schematic  diagram  of  the  microwave  and  spectral  arrangement  used  to 
obtain  radial  profile  data  from  a  low  density  argon  plasma. 
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Figure  2.  Photograph  of  the  R2E-1  test  cell  and  apparatus. 
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Figure  3.  Schematic  sketch  of  the  R2E  1  test  cell. 

2  1  TEST  CELL 

The  viewing  section  ot'  the  lest  eel!  was  constructed  from  four  pieces  of  1  in. -thick 
optically  Hat  quart/.  9  in.  wide  and  10  in.  high  The  l|-in. -square  quart/  test  section  was 
attached  to  a  12.75  in.-diam  section  of  pipe.  The  lower  section  ot  pipe  (approximately 
36  m  long)  housed  the  plasma  arcjet.  The  upper  section  of  pipe  was  connected  to  vacuum 
pumps.  Test  cell  pressures  as  low  as  0.0!  torr  with  no  flow  and  0.6  torr  with  arc  jet 
and  shroud  flow'  could  be  maintained.  A  schematic  diagram  of  the  test  cell  and  arc  jet 
installation  is  shown  in  Fig.  3. 
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2.2  ARC  GENERATOR  AND  PLASMA 


The  plasma  generator  used  in  these  experiments  is  a  Gerdian-type  configuration  and 
is  fully  described  in  Ref.  4.  The  arc  jet  is  mounted  on  a  hydraulic  jack  so  that  various 
axial  positions  of  the  plasma  can  be  viewed  in  the  test  section.  The  arc  jet  (Fig.  4)  is 
equipped  with  a  circular  shroud  made  from  1  /4-in.-diani  steel  tubing  with  several  small 
holes  drilled  around  the  perimeter  to  permit  the  inclusion  of  a  stream  of  argon  around 
the  argon  plasma.  The  shroud  of  argon  makes  the  plasma  more  stable  and  also  reduces 
the  mixing  of  the  plasma  with  the  ambient  air  in  the  test  section.  Typical  operating 
conditions  are  as  follows,  arc  voltage,  25  v;  arc  current,  168  amp;  chamber  pressure,  5 
in.  Hg;  argon  flow  through  arc  jet,  0.32  gm/scc,  argon  flow  through  shroud,  0.28  gm/sec, 
test  cell  pressure,  0.6  torr. 


The  argon  gas  enters  the  arc  chamber  and  expands  from  a  l/4-in.-diam  sonic  nozzle 
into  the  low  pressure  test  cell  after  passing  through  the  arc.  (The  degree  ot  ionization 
of  an  argon  gas  under  the  operating  conditions  outlined  above  ranges  from  1  to  5  percent 
(Ref.  4)).  A  schematic  diagram  of  a  typical  supersonic  argon  plasma  is  shown  in  Fig. 
5.  fhe  supersonic  expansion  is  bounded  by  a  barrel  shock  up  to  the  Mach  disc,  and  the 
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flow  is  subsonic  just  downstream  of  the  Mach  disc.  The  properties  of  the  flow  are 
considerably  different  in  the  region  downstream  of  the  normal  shock  and  t he  region 
upstream  of  the  normal  shock  (Kef.  5).  Most  of  the  microwave  and  spectral  data  were 
obtained  just  upstream  (Station  1)  and  just  downstream  (Station  2)  of  the  Mach  disc. 
(Only  the  data  taken  at  Station  2  are  reported  herein  because  of  the  difficulty  in  calculating 
the  collision  frequencies  in  the  supersonic  region  of  the  plasma  jet.  Also,  the  barrel  shock 
in  the  supersonic  regions  requires  extensive  care  and  preparation  to  ensure  that  reliable 
data  Irom  the  Abel  inversion  are  obtained.) 


Shroud 


Figure  5.  Schematic  diagram  of  an  argon  plasma  free  jet  expanding 
into  a  low  pressure  test  cell. 

2.3  SPECTROMETRIC  ARRANGEMENT  FOR  INTENSITY  MEASUREMENTS 

A  I  in  Jarrell-Ash®  grating  spectrometer  with  a  C/emy-Turner  mount  was  used  for 
the  intensity  measurements.  The  spectrometer  is  equipped  with  curved  slits  and  1180 
groove/mm  grating  bla/.cd  for  maximum  reflection  at  7500A.  A  40 -/a  slit  was  used  which 
gives  an  equivalent  slit  width  of  0.24 A  and  a  resolution  of  0.25 A  An  F-8 ,  35-mm  focal 
length  achromatic  lens  was  mounted  at  the  entrance  slit  of  the  spectrometer  and  focused 
on  the  plasma  centerline 

A  pair  of  plane  front  surface  mirrors,  Ml  and  M2,  were  located  in  the  optical  path 
as  shown  in  Fig.  I .  Mirror  M 1  could  be  rotated  through  a  5-deg  angle  so  that  a  radial 
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scan  of  the  plasma  could  be  obtained.  The  mirror  was  remotely  operated  by  a  synchronous 
reversible  motor.  The  angular  position  was  monitored  with  a  1000-12  linear  potentiometer 
which  served  as  the  x-input  for  the  x-y  plotter  used  to  record  the  intensity  data.  The 
optical  path  from  mirror  Ml  to  plasma  centerline  was  24  in.,  which  provided  a  total 
radiai  scan  of  6  in.  The  small  angle  of  deflection  coupled  with  the  24-in,  focal  length 
very  closely  approximates  a  straight  line  through  the  center  of  the  plasma.  The  measured 
field  of  view  at  the  focal  point  was  approximately  1x2  mm. 

An  RCA  IP28  photomultiplier  was  used  as  the  detector.  The  system  was  calibrated 
for  absolute  intensity  by  placing  a  secondary  standard  tungsten  ribbon  filament  lamp  at 
the  focal  point  of  the  optics  and  recording  the  known  intensity  over  the  range  of  wavelength 
of  interest. 

2.4  35-GHz  PLASMA  ATTENUATION  SYSTEM 

A  schematic  diagram  of  the  microwave  attenuation  system  and  the  hydraulically 
operated  table  for  scanning  the  plasma  is  given  in  Fig.  1 .  A  forced-air-cooled  reflex  klystron 
operating  at  400  mw  at  a  frequency  of  34.91  GHz  provided  the  microwave  power  for 
the  experiments.  The  microwave  signal,  after  leaving  the  klystron,  passes  through  a  ferrite 
isolator  to  reduce  the  strength  of  any  signals  which  might  be  reflected  back  into  the 
klystron.  A  0-  to  -20-db  attenuator  was  added  so  that  the  power  to  the  antenna  could 
be  adjusted  as  required.  (Generally  this  attenuator  was  set  at  the  zero  position).  The  power 
output  was  monitored  by  including  a  20-db  directional  coupler  and  monitoring  the  output 
signal  (SI)  with  a  crystral  detector.  (If  the  klystron  and  power  supply  were  turned  on 
2  hrs  prior  to  testing,  it  was  found  that  the  klystron  maintained  a  steady  output  with 
no  drift.  Therefore,  a  2-m  warmup  period  was  always  observed  prior  to  testing).  Following 
the  20-db  coupler  is  a  0-  to  55-db  precision  attenuator  which  was  used  for  calibration 
of  the  microwave  attenuation  system.  The  calibration  extended  from  0  to  45  db.  Another 
directional  coupler  was  placed  after  the  attenuator  and  was  used  to  monitor  the  reflected 
signal  (S2).  (The  reflected  signal  mentioned  here  represents  the  reflected  signals  from  all 
points  beyond  the  coupler  (e.g.,  the  transition  guide,  antenna,  lens,  test  cell,  etc.),  and 
it  serves  primarily  as  an  aid  in  adjusting  the  E/H  tuner  for  impedance-matching  purposes, 
and  not  as  measure  of  the  reflected  signal  from  the  plasma).  An  E/H  tuner  was  inserted 
prior  to  entering  the  antenna  for  impedance-matching  purposes.  The  signal  was  then  fed 
to  a  dielectric  lens-focused  antenna,  then  through  the  test  cell  and  into  the  receiving  antenna 
from  which  the  attenuated  signal  (S3)  was  monitored. 

The  focusing  antenna  consisted  of  two  confocal  conical  horns  (7.032  in.  long  with 
an  aperture  of  6  in.,  phase  corrected  with  dielectric  lens,  Rexolite®  1422  (Refs.  6  and 
7)).  The  horn  antennas  were  originally  designed  for  X-band  operation;  therefore,  a  transition 
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section  from  K-  to  X-band  was  required.  The  antennas  used  in  these  experiments  were 
not  designed  specifically  for  optimum  conditions  (i.e.,  the  flare  angle  of  the  antenna  for 
either  the  E  or  H  plane  was  not  chosen  for  optimum  operating  conditions  (see  Refs. 
8  and  9).  Also,  the  energy  profile  emerging  from  the  waveguide  is  not  uniform.  These 
factors  result  in  an  uneven  illumination  of  the  fust  lens,  and  the  net  effect  is  that  the 
energy  of  the  beam  is  concentrated  in  a  narrow  beam  of  approximately  1-cm  diameter 
with  a  fairly  long  focal  depth.  Thus,  the  system  was  insensitive  to  minor  variations  in 
focal  length.  This  feature  actually  simplified  the  laboratory  operation  and  the  analysis 
because  the  main  body  of  the  microwave  energy  was  concentrated  along  the  microwave 
antenna  centerline  in  a  narrow  beam  of  plane  waves  as  they  passed  through  the  plasma. 

3.0  MICROWAVE  ATTENUATION  CALIBRATION 

When  reflections  from  plasma  boundaries  and  internal  refractions  can  be  neglected, 
then  the  plasma  can  be  treated  as  an  insertion  loss,  and  the  signal  at  the  receiving  antenna 
can  be  calibrated  by  inserting  calibrated  absorbtion  devices  into  the  signal  path.  (Note: 
reflections  can  be  neglected  if  the  plasma  properties  vary  slowly  with  respect  to  the 
wavelength  of  the  incoming  microwave  energy;  see  Ref.  10.)  The  microwave  system  was 
calibrated  in  these  experiments  by  inserting  a  precision  attenuator  in  the  system  (see  Fig. 
1)  and  recording  the  output  signal  (S3)  as  the  attenuator  was  varied  from  0  to  45  db. 
The  large  range  of  attenuation  represents  a  power  variation  of  approximately  five  orders 
of  magnitude,  and  therefore  requires  a  great  deal  of  care  at  the  higher  attenuation  values 
in  order  to  maintain  an  accurate  recording  of  these  calibration  points.  A  high  gain  amplifier 
was  added  at  the  receiving  detector  to  amplify  the  signal  (S3)  at  the  higher  values  of 
attenuation.  The  addition  of  the  amplifier  required  that  the  calibration  be  performed  at 
several  (typically  five)  known  gain  settings  and  that  the  plasma  attenuation  data  be  obtained 
at  the  same  gain  settings. 

The  operating  procedures  for  calibration  were  to  set  the  gain  of  the  amplifier  and 
then  to  set  the  precision  attenuator  to  a  known  value.  The  attenuated  signal  (S3)  and 
position  signal  (S6)  were  recorded  on  the  x-y  plotter  as  the  microwave  beam  was  scanned 
across  the  test  section  with  no  plasma.  The  attenuated  signal  was  recorded  as  the  y-input 
on  an  x-y  plotter,  and  the  signal  which  monitored  position  of  the  table  (S6)  was  used 
as  the  x-input.  The  gain  was  held  constant  and  this  procedure  repeated  for  a  number 
of  attenuation  settings  on  the  precision  attenuator.  The  gain  was  then  increased  to  a  new 
level  and  the  above  procedure  repeated  until  calibration  data  were  obtained  over  the  desired 
range  (typically  0  to  45  db).  The  precision  attenuator  was  then  reset  to  0  db  in  preparation 
for  plasma  attenuation  measurements. 
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4.0  ABEL  INVERSION  TECHNIQUE 

The  measured  intensity  of  an  emission  line  in  an  optically  thin  plasma  represents 
the  integrated  sum  of  the  emissions  of  the  emitting  species  along  the  optical  path  in  the 
plasma.  Similarly,  the  microwave  signal  at  the  receiver  end  represents  the  transmitted  signal 
minus  the  integrated  loss  of  the  transmitted  signal  as  it  passes  through  the  plasma.  If 
the  property  being  observed  (i.e..  emission  coefficient,  absorption  coefficient,  etc.)  is 
constant  along  the  observation  path,  the  integral  is  simply  the  product  of  the  coefficient 
times  the  path  length.  However,  when  the  properties  are  not  uniform  along  the  path  of 
observation,  an  inversion  scheme  is  required  to  obtain  local  values  of  the  properties. 

The  Abel  transformation  (Ref.  11)  represents  a  mathematical  procedure  by  which 
the  integrated  signal  from  a  cylindrically  symmetric  source  can  be  transformed  into  local 
properties.  The  main  features  of  this  technique  are  illustrated  in  Fig.  6.  If  it  can  be  assumed 


y 


Figure  6.  Schematic  diagram  of  a  cylindrically  symmetric  system  illustrating 
the  coordinate  notation  used  in  the  Abel  inversion. 

that  the  measured  parameter  G(y)  represents  the  integrated  sum  of  the  local  property 
f(r),  then  the  following  analysis  may  be  made. 

\.>  X  r 

G(v)  -  I  f(r)dx  -  2  |  f(r>dx  (!) 

--v ,  <' 
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where 


x 


and 


/o2  2 

x,  -  vK  -v 


The  transformation  to  cylindrical  coordinates  yields 


H  ((r)rdr 

«y>  -  2!  -r t=T 

7  \i  T  —  V 


(2) 


(3) 


Thus  under  the  conditions  outlined  above,  the  local  plasma  parameter,  f(r),  can  be 
obtained  using  the  Abel  transformation  as  shown  below. 


]  r  ^  (rtO/Jvbly 

n  ?  f~2 - 2 

V  >'  -  r 


C4) 


Many  techniques  have  been  developed  to  perform  the  inversion  calculation  of  Eq. 
(4).  The  choice  of  technique  depends  primarily  upon  the  investigator's  preference,  and 
to  some  extent  upon  the  quality  of  the  data.  The  procedure  for  evaluating  Eq.  (4)  in 
this  investigation  was  to  represent  the  measured  parameter  G(y)  (either  line  intensity  or 
attenuation)  by  a  series  of  least-square-fit,  even-function,  sixth-order  polynominals  and 
to  perform  the  differentiation  and  integration  analytically  using  numerical  methods.  The 
validity  of  the  curve  fit  requires  data  fluctuations  to  be  random.  Consequently,  it  is 
necessary  that  several  data  points  at  each  radial  position  be  input  to  the  computer  program 
in  order  to  obtain  an  accurate  least-square  fit  to  these  random  data  points.  However, 
when  the  fluctuations  in  the  data  are  small  (i.e.,  when  the  data  can  be  characterized  by 
a  single  curve),  then  the  curve  fit  must  be  exact.  Otherwise,  large  fluctuations  in  the 
inverted  data  can  be  expected.  Since  the  data  used  in  this  report  were  of  the  latter  category 
(i.e..  exhibiting  very  small  fluctuations),  it  was  necessary  to  demand  a  close  curve  fit  (±0.5 
percent)  to  the  data  points.  It  was  possible  to  obtain  the  curve  lit  within  ±0.5-percent 
tolerances  only  from  the  centerline  out  to  approximately  1  cm  and  within  ± 5-percent 
tolerance  on  out  to  10  or  11  cm. 


As  can  be  seen  from  Eq.  (4)  the  coefficient  f(r)  is  functionally  related  to  the  spatial 
derivative  of  the  measured  parameter  G(y),  and  it  is  therefore  possible  to  produce  major 
fluctuations  in  f(r)  with  only  minor  fluctuations  in  G(y),  which  is  a  main  weakness  of 
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the  inversion  technique.  In  fact,  it  requires  care  and  skill  in  order  to  properly  apply  any 
inversion  scheme  to  real  cases  because  of  the  sensitivity  of  the  inverted  results  to  the 
quality  of  the  measured  data. 

5.0  ELECTRON  TEMPERATURE  AND  NUMBER  DENSITY 
FROM  SPECTRAL  MEASUREMENTS 


The  electron  temperature  and  density  of  an  electrically  neutral  plasma  for  which 
the  upper  energy  states  of  the  excited  argon  are  in  a  Boltzmann  distribution  are  related 
according  to  Saha's  equation  (Ref.  12). 


<7 


kTe(r>y/2N 

.2  /  n  -jeIcc 


Z 


exp  [  -  ( Ej  ~  En)  T  Tft(r)] 


(5) 


For  a  given  atom  there  are  two  unknowns  in  this  equation;  Nn  (the  number  of  atoms 
which  are  in  the  energy  state  En)  and  Te  (the  electron  temperature).  Thus  two  parameters 
must  be  measured  to  obtain  the  electron  number  density1 .  The  measurement  of  the 
intensity  of  two  emission  lines  is  sufficient  to  obtain  these  parameters.  That  is,  the  power 
radiated  per  unit  volume  per  unit  solid  angle  at  a  given  wavelength  is 

Sin  ”  l-X  (6) 


where  enm  is  often  called  the  emission  coefficient,  X  is  the  wavelength,  and  Anm  is  the 
transition  probability  associated  with  the  transition  to  En  (Ref.  II).  The  intensity  of  the 
radiation  at  the  spectrometer  is 


in„/y>  =  2/ 


nm 


Use  of  the  Abel  transform  (see  Eqs.  (4)  and  (6))  yields 


4A  fr 

hc  ''nm  H  ^  2  _  2 


(7) 


(8) 


*  If  the  energy  states  were  in  equilibrium  down  to  the  ground  state,  it  would  be  possible  to  calculate  Nn  from 
pressure  and  election  temperature  data.  For  low  density  Argon  plasmas  this  is  seldom  true. 
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The  ratio  of  the  emission  coefficients  of  two  argon  lines  labeled  1  and  2  can  be 
used  to  obtain  the  electron  temperature  (Ref.  11). 


where  E  is  the  energy  of  the  upper  state,  k  is  Boltzmann's  constant,  g  is  the  statistical 
weight,  A  is  the  transition  probability,  and  X  is  the  wavelength.  The  emission  coefficients 
of  the  radiation  for  the  two  different  emission  lines  subscripted  1  and  2  are 


and 


,<r>  =  i  / 


(dl ]/ dy)dy 

VI  2 

y  -  r 


(10) 
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6.0  SPECTRAL  AND  MICROWAVE  MEASUREMENTS 


(11) 


6.1  SPECTRAL  CONSIDERATIONS 

Originally  it  was  planned  to  use  three  or  four  argon  emission  lines  in  order  to  obtain 
an  average  value  of  temperature  and  electron  density  at  each  radial  position  in  the  plasma. 
In  practice,  however,  this  proved  to  be  difficult  because  of  the  low  intensity  of  some 
lines  of  interest  which  originated  from  states  near  the  continuum,  and  also  because  it 
was  necessary  to  keep  the  measuring  time  to  a  minimum  in  order  to  reduce  the  possible 
variations  in  the  data  which  might  result  from  plasma  and  instrumentation  drift.  The 
spectral  response  of  the  photomultiplier  also  constrains  the  selection  of  useful  emission 
lines. 


The  41 58-,  6032-,  and  6965-A  lines  were  chosen.  It  was  subsequently  discovered, 
however,  that  the  assumption  of  optically  thin  plasma  and/or  Boltzmann  distribution  was 
not  applicable  for  the  6965-A  line.  The  characteristic  flatness  of  the  intensity  profile  in 
the  data  obtained  with  this  line  indicated  that  self-absorption  might  be  occurring.  The 
low  values  of  electron  temperature  consistently  obtained  when  the  data  from  the  6965-A 
line  were  used  for  temperature  calculations  generated  doubt  as  to  whether  the  upper  energy 
state  associated  with  this  line  was  in  a  Boltzmann  equilibrium  distribution  with  the  other 
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upper  energy  states,  Thus  the  6032-  ami  4158-A  lines  were  used  in  this  report  to  obtain 
electron  temperature  and  electron  density  profiles. 

6.2  SPECTRAL  MEASUREMENTS 

The  spectrometer  was  carefully  adjusted  so  that  the  peak  intensity  of  each  line  of 
interest  was  established,  and  then  the  plasma  was  scanned  by  rotating  the  rotating  mirror 
Ml  (see  Fig.  1)  and  recording  the  photomultiplier  output  (S4)  on  an  x-y  plotter.  The 
plasma  was  scanned  several  times  for  each  line  being  observed  to  ensure  repeatability. 
This  procedure  was  repeated  for  each  emission  line  of  interest.  The  measured  intensity 
data  obtained  were  used  in  the  Abel  inversion  program,  for  which  the  emission  coefficient 
for  each  line,  the  electron  temperature,  and  electron  number  density  were  calculated. 

The  measured  intensity  profiles  of  the  4158-  and  6032-A  Lines  taken  at  Station  2 
(approximately  2  in.  downstream  of  the  Mach  disk)  are  shown  in  Fig.  7.  The  error 
boundaries  (indicated  by  the  dashed  lines  in  Fig.  7)  were  obtained  from  the  error  bounds 
quoted  for  the  calibration  curve  which  was  used  to  convert  the  recorded  photomultiplier 
voltage  to  radiant  intensity  units.  Instrumentation  errors  were  neglected  because  the  same 
instruments  were  used  under  identical  circumstances  for  intensity  calibration  and  for  plasma 
spectral  measurements.  The  inverted  emission  coefficients,  with  the  indicated  error  bounds 
for  these  two  coefficients,  are  shown  in  Fig.  8.  The  error  bounds  shown  in  Fig.  8  for 


Cartesian  Coordinate  y,  cm 

Figure  7.  Measured  intensity  profiles  for  the  4158-  and  6032-A  argon 
emission  lines. 
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Figure  8.  Emission  coefficient  profiles  for  the  4158-  and  6032-A  argon 
emission  tines. 


the  emission  coefficients  are  the  same  as  those  for  the  intensity  curve;  that  is,  it  is  assumed 
that  the  inversion  process  does  not  generate  any  additional  error,  and  therefore,  that  the 
error  bounds  should  be  the  same  for  both  the  intensity  and  the  emission  coefficients. 
The  electron  temperature  and  electron  number  density  using  4158-  and  6032-A  emission 
data  lines  are  shown  in  Figs.  9  and  10,  respectively.  The  principal  factors  which  contribute 
to  the  uncertainty  in  the  electron  temperature  are  the  transition  probability  uncertainties 
(  +  5  percent  for  the  4158-A  transition  and  ±5  percent  for  the  6032-A  transition  (Ref. 
13))  and  the  uncertainty  associated  with  the  emission  coefficients  (±4  percent).  (The 
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Radial  Coordinate,  cm 


Figure  9.  Plasma  electron  temperature  profile  from  line  ratio  calculations. 


± 4-percent  uncertainty  in  the  emission  coefficients  is  based  upon  the  uncertainty  associated 
with  the  secondary  tungsten  filament  used  for  absolute  intensity  calibrations  (±3  percent) 
and  reproduction  uncertainty  (±1  percent)).  The  errors  combine  to  give  an  uncertainty 
in  the  electron  temperature  obtained  for  line  ratio  calculations  of  approximately  ±  1 8 
percent.  The  uncertainty  in  the  electron  number  density  obtained  from  absolute  line 
intensity  measurements  can  be  traced  to  the  uncertainty  in  the  transition  probability  (±5 
percent),  emission  coefficients  (±4  percent),  and  electron  temperature  (±18  percent).  The 
net  effect  of  these  errors  is  to  produce  an  uncertainty  of  ±36  percent  in  the  electron 
density  obtained  from  spectroscopic  data. 
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Figure  10.  Plasma  electron  density  profits  using  absolute  intensity 
measurements  of  the  4158-  and  6032-A  argon  emission 
Ifties. 
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6.3  MICROWAVE  ATTENUATION  MEASUREMENTS 
6.3.1  Measuring  Procedures 

The  test  cell  was  evacuated  and  the  microwave  system  was  calibrated  as  outlined 
in  Section  3.  Next  the  plasma  flow  was  initiated,  and  the  plume  was  positioned  in  the 
test  cell  so  that  the  desired  axial  position  along  the  plasma  could  be  viewed.  The  microwave 
attenuation  data  were  then  recorded  on  an  x-y  plotter  as  the  microwave  beam  was  scanned 
radially  across  the  plasma.  This  procedure  was  repeated  tor  each  value  of  amplifier  gain 
used  during  the  system  calibration.  The  attenuation  profile  was  obtained  by  superimposing 
the  calibration  curves  onto  the  measured  data.  A  representative  attenuation  profile  obtained 
by  this  procedure  is  shown  in  Fig.  11.  The  overlapping  data  points  observed  result  from 
the  overlap  obtained  at  adjoining  gain  positions  used  in  calibration  and  data  acquisition. 


Cartesian  Coordinate  x,  cm 


Figure  11.  Plasma  microwave  attenuation  data  2  in.  downstream  of 
normal  shock. 

The  uncertainty  in  the  attenuation  data  can  be  traced  to  the  uncertainty  in  the  precision 
attenuator  (Fig.  1)  which  was  used  for  calibration  purposes  and  to  the  random  error 
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generated  in  the  data  acquisition  process.  The  total  error  associated  with  the  attenuation 
data  is  estimated  to  be  approximately  ±5  percent. 

6.3.2  Attenuation  Coefficient 

The  attenuation  profile  shown  in  Fig.  1 1  represents  the  integrated  attenuation  of 
the  microwave  signal  as  it  passes  through  the  plasma.  Many  authors  (e.g.  see  Refs.  14, 
15.  and  16)  have  used  a  rather  simple  relationship  between  the  measured  integrated 
attenuation  and  the  attenuation  coefficient  (see  Appendix  A),  and  in  these  cases  the  process 
of  obtaining  the  inverted  data  becomes  a  straightforward  procedure.  That  is,  ordinary 
Abel  inversion  techniques  are  applicable.  Thus,  when  it  is  assumed  that  the  microwave 
power  Pout  exiting  the  absorbing  medium  is  given  by 

pout  =  pin  exP  t-  /«<*>']  (12) 

(where  PIN.  is  the  input  microwave  power,  and  the  line  integral  of  the  attenuation 
coefficient  is  taken  over  the  ray  path-see  Ref.  16),  then  the  attenuation  in  units  of  decibels 
becomes 


A( k )  =  10  log10\"7r~J  =  -  4.343  fady 
\  IN/  i 

For  a  cylindrically  symmetric  plasma,  Eq.  (13)  becomes 


(13) 


A( x)  =  -  (2)(4.343) />  a(r)dy  (14) 

O 

It  is  noted  that  Eq.  (14)  is  of  the  standard  Abel  form,  and  therefore  the  attenuation 
coefficient  can  be  obtained  quite  easily  (see  Section  4).  That  is, 


a(r) 


(dA/dx)djc 

'  yfJZ~ 


(15) 


The  numerical  procedure  for  evaluating  this  expression  is  outlined  in  Section  5.  The  radial 
profile  of  the  attenuation  coefficient  for  the  measured  integrated  attenuation  shown  in 
Fig.  II  for  the  conditions  outlined  above  is  presented  in  Fig.  12. 
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Figure  12.  Microwave  attenuation  coefficient  versus  radius. 


6.3.3  Collision  Frequency 

In  order  to  calculate  the  electron  density  from  the  attenuation  coefficient  it  is  also 
necessary  to  either  measure  or  calculate  the  electron  collision  frequency.  For  a  pure  argon 
plasma  there  is  sufficient  data  extant  (Ref.  17)  so  that  reliable  equations  can  be  used 
to  calculate  the  collision  frequency  (Ref.  7).  Only  two  types  of  encounters  need  be 
considered:  (1)  encounters  between  electrons  and  neutral  argon  atoms,  and  (2)  encounters 
between  electrons  and  argon  ions.  (Electron-electron  encounters  do  not  change  the 
transport  of  charge,  and  therefore  may  be  ignored).  That  is,  the  collision  frequency  v  can 
be  written  as  the  sum  of  the  electron-neutral  collision  frequency  ven  and  the  electron-ion 
collision  frequency  (Ref.  10).  Thus, 
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V 


(16) 


where 


1/  =  2.6  x  lO4  Q  n  T 

en  en  g  e 


1/2 


(17) 


and 


vei  =  3.63  x  10*6(noTe-3/2)fnA 
A=  1.23  x  104  Te3/2  ne‘I/2  (18) 

Values  for  the  collision  cross  section,  Qen>  as  a  function  of  electron  temperature,  Te, 
were  obtained  from  Ref.  17.  (The  electron  temperature  was  determined  from  the  spectral 
data).  The  gas  density,  ng,  at  Station  2  was  obtained  by  assuming  that  the  static  pressure 
in  the  stream  is  equal  to  the  test  cell  pressure  and  assuming  that  the  gas  temperature 
is  approximately  equal  to  the  total  temperature  (Ref.  5). 

6.3.4  Electron  Density 

The  relationship  between  the  electron  density,  ne,  and  the  attenuation  coefficient, 
a,  is  as  follows  (see  Appendix  A): 

(19> 

where  c  is  the  speed  of  light,  o>  is  the  microwave  frequency,  me  is  the  electron  mass, 
and 
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(20) 


(21) 


The  electron  number  density,  n.,  may  be  calculated  directly  from  Eq.  (19)  using 
the  relationships  established  between  ne,  PR ,  and  Pj  in  Eqs.  (20)  and  (21).  That  is, 
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n  - 

i 


(22) 


In  the  limit  v1  «  cop2  «  to2  (where  coP  =  (nce2/mee0)1'2  =  plasma  frequency),  an 
approximate  expression  for  a  can  be  obtained  from  Eq.  (19)  (Ref.  14).  That  is, 


a 


•t.62  x 


10'S 


n  i' 

V 

~2  ~2 
fit  t  J' 


Solving  for  nc. 


n 


2.165  x 


(23) 


(24) 


Since  the  collision  frequency  is  also  a  function  of  the  electron  density,  it  is  necessary 
to  use  an  iterative  scheme  to  obtain  nc  versus  a  (see  Section  6.3.3).  Equations  (22)  and 
(24)  are  plotted  in  Fig.  13  for  various  values  of  v.  In  this  figure  it  can  be  seen  that 
the  approximate  expression  for  ne  (see  Eq.  (24))  asymptotically  approaches  the  exact 
expression  (see  Eq.  (22))  when  cu  »  top  for  values  of  v  <  co.  However,  for  <o  <  cop 
the  approximate  expression  differs  greatly  from  the  exact  expression. 


Figure  13.  Electron  density  calculations  versus  attenuation  coefficient 
for  various  values  of  collision  frequency. 
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7.0  DISCUSSION  OF  RESULTS 


The  radial  profiles  of  electron  density  obtained  from  the  spectroscopic  data  (Fig. 
10)  and  the  profiles  obtained  from  the  microwave  data  in  conjunction  with  Eqs.  (22) 
and  (24)  are  compared  in  Fig.  14.  The  error  bounds  for  the  electron  densities  obtained 


0  1  2  3  4  5 


Radial  Coordinate,  cm 


Figure  14.  Comparison  of  plasma  electron  density  profiles  using 
microwave  attenuation  and  spectral  data. 
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from  microwave  attenuation  measurement  result  from  uncertainties  in  the  value  of  the 
absorbtion  coefficient  and  the  electron  temperature.  The  error  bounds  for  the  electron 
density  are  ±33  percent  using  Eq.  (22)  and  ±26  percent  using  Eq.  (24).  The  plasma  electron 
density  profiles  obtained  from  the  microwave  measurements  (either  Eq.  (22)  or  Eq.  (24)) 
do  not  agree  with  the  electron  density  obtained  from  spectral  measurements.  The  difference 
between  the  profiles  obtained  from  microwave  measurements  (see  Eqs.  (22)  and  (24)) 
can  be  explained  by  noting  that  for  the  plasma  conditions  at  Station  2  the  limits  specified 
in  the  simplified  expression  for  electron  density  (Eq.  (24))  are  not  applicable. 

The  values  of  the  electron  temperature,  Te,  total  collision  frequency,  v ,  electron-ion 
collision  frequency,  vcl.  angular  plasma  frequency,  toP,  electron  density,  ne  (Eq.  22),  and 
the  values  of  electron  density  obtained  from  spectral  measurements  at  various  radial 
positions,  r,  are  tabulated  for  the  argon  plasma  at  Station  2  (Fig.  5)  in  Table  1 .  Examination 
of  this  table  shows  that  the  limits  required  for  Eqs.  (23)  and  (24)  to  be  applicable  are 
not  satisfied. 

Table  1.  Electron  Temperature,  Total  Electron  Collision  Frequency,  Electron-Ion 
Collision  Frequency,  Plasma  Frequency,  and  Electron  Densities  at  Various 
Radial  Positions  Approximately  2  in.  Downstream  of  the  Mach  Disk 


r. 

Radial 

Position, 

Ci:i 

Te, 

Election 

Temperature. 

K 

V, 

Total  Collision 
Frequency, 

«c"*  a  10® 

"«• 

Electron-Ion 
Frequency. 
sec-1  a  10® 

r 

Plasma  Frequency, 
nd/sec  s  I01 1 

nc. 

Eq  (22). 
cm"3  x  I013 

nc. 

Spectral. 
enr3  x  10' 3 

oo 

73i2 

7.14 

707 

2  27 

l  62 

9  35 

05 

7175 

6  95 

6  87 

2  24 

l  58 

4. IS 

1  0 

7600 

6.52 

6.44 

221 

1.54 

8  70 

1.5 

7880 

6  19 

6  11 

2  21 

1  53 

8  05 

20 

8225 

5  84 

5.75 

2  20 

1  S2 

6  92 

25 

8375 

5  70 

5.61 

2.20 

1  52 

5.20 

3.0 

7480 

6.57 

6.49 

2  20 

1  52 

2  55 

35 

6150 

7.89 

7.83 

2  20 

1.52 

0.89 

4.0 

5150 

9  79 

9.74 

2  20 

]  .5 1 

0.74 

4.5 

4850 

1230 

12.23 

2  19 

1 

0.4 1 

5.0 

3820 

15.20 

15.14 

2  17 

1.47 

0.19 

5.5 

3200 

18  50 

18.49 

212 

]  41 

0.10 

The  difference  between  the  spectral  results  and  the  microwave  results  (particularly 
those  obtained  using  the  exact  expression  given  by  Eq.  (22))  are  more  difficult  to  explain. 
It  is  possible  that  the  results  obtained  from  spectral  data  are  not  correct.  However,  a 
more  likely  explanation  of  the  difference  is  that  the  microwave  equations  (Eqs.  (22)  and 
(24))  which  are  used  to  relate  the  attenuation  coefficient  to  the  integrated  total  attenuation 
are  not  applicable  to  inhomogeneous  media.  First  of  all,  Eq.  (13)  is  based  upon  a  plane 
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wave  passing  through  a  uniform  conductive  media  (see  Appendix  A)  and  not  upon  the 
case  in  which  V<r=£o  (i.e.,  the  propagation  coefficient  in  the  wave  equation  is  not  constant). 

Since  the  microwave  beam  cross  section  (approximately  1-cm  diameter)  is  much  larger 
than  the  spectral  cross  section  (1x2  mm),  it  is  reasonable  to  expect  that  some  differences 
in  the  density  profiles  obtained  from  microwave  and  spectral  data  would  exist.  That  is, 
the  larger  aperture  associated  with  the  microwave  system  is  more  effective  at  averaging 
out  sharp  changes  in  the  density  profile.  Although  the  averaging  effect  could  account 

for  some  of  the  discrepancies  between  the  density  profile  obtained  from  spectral 

measurements  and  the  density  profile  obtained  from  microwave  measurements  (Eq.  (22)) 
(particularly  toward  the  outer  boundaries  of  the  plasma),  averaging  cannot  account  for 
the  large  difference  (approximately  a  factor  of  five)  at  the  plasma  centerline. 

Other  sources  of  microwave  error  can  arise  from  reflections  which  occur  at  sharp 

boundaries  and  from  signal  scattering  in  the  plasma.  A  cylindrical  plasma  in  conjunction 

with  a  confocal  microwave  system  makes  it  extremely  difficult  to  measure  the  reflected 
signal  because  of  the  geometry  of  the  relecting  angles.  However,  unless  the  plasma 
boundaries  are  very  sharp  (i.e.,  unless  electron  density  changes  rapidly  in  distances  of 
less  than  a  wavelength),  boundary  reflections  should  be  minimal  (Ref.  10). 

The  express  purpose  of  this  project  was  to  make  simultaneous  measurements  of  the 
plasma  electron  density  profile  using  spectral  data  and  microwave  data  with  the  intent 
of  checking  the  microwave  calculations  against  the  calculations  obtained  from  spectral 
data.  The  primary  purpose  was  not  to  resolve  any  existing  differences  but  rather  to  provide 
reliable  experimental  data  from  which  further  studies  might  be  obtained. 

8.0  SUMMARY  OF  RESULTS 

The  special  arrangement  of  the  microwave  and  spectral  instrumentation  along  with 
the  test  cell  configuration  offered  an  ideal  opportunity  to  obtain  reliable  spectral  and 
microwave  attenuation  data.  For  the  most  part,  the  integrated  intensity  profiles  and  the 
microwave  attenuation  profiles  were  symmetrical  about  the  axis  of  the  plasma;  therefore, 
the  inversion  of  the  data  could  be  obtained  by  straightforward  procedures. 

The  calculation  of  the  electron  density  profiles  from  these  data  (spectral  versus 
microwave)  did  not  agree.  The  large  error  bounds  for  the  calculations  of  electron  number 
density  indicate  the  problem  of  attempting  to  make  a  clear-cut  comparison  between 
electron  density  calculations  based  on  microwave  attenuation  and  spectral  line  emission 
data.  The  large  error  bounds  for  the  electron  density  calculations  are  based  upon  the 
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uncertainty  in  the  basic  physical  constants  associated  with  transitions  between  excited 
states  for  argon.  The  uncertainty  in  the  collision  frequency  is  the  principle  factor  which 
contributes  to  the  uncertainty  in  the  microwave  electron  density  calculations.  The  collision 
frequency  uncertainty  can  be  traced  to  the  uncertainty  in  the  electron  temperature,  and 
therefore,  the  calculations  of  electron  density  from  microwave  attenuation  calculations, 
in  this  instance,  are  also  subject  to  the  uncertainty  attributed  to  basic  physical  constants 
of  the  argon  emission  lines.  The  situation  becomes  even  worse  when  the  plasma  is  composed 
of  gases  more  complex  than  an  argon  plasma.  Other  methods  of  obtaining  the  collision 
frequency  for  more  accurate  electron  density  measurements  by  microwave  attenuation 
techniques  are  needed. 
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APPENDIX  A 

PLANE  ELECTROMAGNETIC  WAVES  IN  A  HOMOGENEOUS  PLASMA 


Maxwell's  equations  in  an  isotropic  plasma  and  for  time-harmonic  variations  of  the 
form  exp  (-jcjt)  may  be  written 


V  x 

F  = 

jo,  13 

►  0  » 

L1  L  *■  L* 

(A-l) 

V  X 

1  = 

r- 

! 

clr 

1 

..  .  1  > 

e  . 

(A-2) 

- * 

V  • 

K  = 

0 

(A-3) 

V  • 

B  - 

0 

(A-4) 

where  j  =_V-T  and  k  =  (cjfyQe0  +  jcojUoff)1  ^  propagation  coefficient.  The  constitutive 
equation,  D  =  e0E.  13  =  and  Ohm's  law  J  =  crE  were  used  to  obtain  this  particular 
form  of  Maxwell's  equations.  The  free-spaee  values  of  permittivity  e0  =  1/36jt  x  1 0  ^ 
F/m  and  the  free-space  value  of  permeability  =  4?r  x  10' 7  hy/m  were  used  for  the 
plasma.  The  equation  for  the  propagation  of  the  electric  field  through  the  plasma  is 
obtained  by  raking  the  curl  of  Eq.  (A-l)  and  eliminating  B  using  Eq.  (A-2).  That  is, 

v  x  (\v  r!»  -  yrv.  d  _  v2f:  -=  k2F  (A-5) 

— ►  ^ 

Since  V  •  E  =  0  in  a  plasma,  the  wave  equation  for  the  electric  field  propagating  through 
a  plasma  is 


V-K  -  k“F.  =  0 


(A-6) 


Similarly,  the  wave  equation  for  the  magnetic  field  is  obtained  by  taking  the  curl 
of  Eq.  (A-2).  Thus 


V  x  <Y  -  13)  =  S<\  .  B)  -  \“2l1 


(A-7) 


where 


\-  <  kK  -.  \,2<V  <  IM  -  i  x  Vic2  (A-8) 

The  wave  equation  for  the  magnetic  field  is  obtained  from  Eqs.  (A-2),  (A-4).  (A-7).  and 
(A-8).  That  is. 
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V2l1  +  k2B  +  — 


K) 


=  0 


!<4> 


( A-9 ) 


where 


Vic 2  = 


A<7 


Thus.  Eq.  (A-9)  can  be  written  as 


V2B  -  k2B  +  njvc  v  £)  =  o  (A-10) 

It  can  be  seen  from  Eq.  (A-10)  that  when  the  gradient  of  the  conductivity  is  not  zero, 
the  equations  for  the  electric  and  magnetic  fields  are  coupled  and  nonlinear.  A  nonzero 
conductivity  gradient,  of  course,  compounds  the  problem  of  calculating  the  attenuation 
of  the  electromagnetic  field  as  it  travels  through  an  inhomogeneous  plasma  because  it 
becomes  necessary  to  know  not  only  the  value  of  k,  but  also  how  k  varies  from  point 
to  point  in  the  plasma.  However,  before  pursuing  the  problem  of  calculating  the  fields 
and  attenuation  for  the  inhomogeneous  case,  it  is  instructive  to  first  look  at  the 
homogeneous  case. 

It  is  convenient  to  divide  the  complex  propagation  coefficient,  k.  into  real  and 
imaginary  components.  Before  this  is  done,  however,  it  is  necessary  to  note  that  the 
conductivity  is  also  complex.  For  a  Lorentz  plasma  (i.e.,  one  where  the  collision  frequency 
is  independent  of  velocity1),  the  electrical  conductivity  is  of  the  form  (Ref.  14) 


a 


•'  Jf/I 


Thus  the  propagation  coefficient,  k,  is 


k 


to/£oaj) 


i(onrpnY'i  =  P  f  ja 


(A-ll) 


(A- 12) 


where  a  is  the  attenuation  coefficient  and  p  is  the  phase  coefficient.  Solving  for  a  and 
j3  in  Eq.  (A-12)  yields 


a 


-  (1 


-p'iF 


<A-13) 


Ifhe  velocity  here  is  the  avciagc  directed  velocity  for  the  electron  ensemble  and  not  the 


rms  random  velocity. 
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and 


/3  =  R/nT 
V  - c  L 


P,)J 


1  n 


n 


-p>r 


(A- 14) 


where 


P,  =  a  |  W0  and  =  a\Vwti 

The  terminology  used  to  denote  a  and  (3  becomes  evident  when  plane  wave  solutions 
to  the  wave  equations  (Eqs.  (A-6)  and  (A-10))  are  obtained  for  a  homogeneous  plasma. 
The  electric  and  magnetic  fields  for  a  plane  electromagnetic  wave  traveling  in  the  z  direction 
in  a  homogeneous  plasma  (e.g.,  Vo  =  0)  which  are  solutions  to  the  wave  equations  are 

E(z,t)  =  i  jEti  exp \~at  —  j(er>t  ~  j0z)]  (A-15) 


and 


B(z,t) 


exp  t-oz  -  j(uji  -  /3z  -  0)\ 


(A- 16) 


where  E0  is  value  of  the  electric  field  at  z  =  0  and  8  -  tan -l(a//3)  (Note:  only  the  real 
part  of  these  expressions  is  used  in  computing  the  values  of  E  or  B.)  Equations  (A- 1 3) 
and  (A-16)  represent  solutions  to  the  wave  equations  in  which  the  fields  are  attenuated 
by  the  factor  exp(-az)  and  phase  shifted  by  the  factor  exp  (-j0z)  as  the  electromagnetic 
wave  propagates  through  the  plasma  in  the  z  direction. 


The  attenuation  of  a  microwave  signal  as  it  passes  through  a  lossy  medium  in  units 
of  decibels  (db)  is  defined  as 


A(dh)  - 


(A-l  7) 


where  Poul  is  the  electromagnetic  power  emerging  from  the  lossy  medium  and  Pin  is 
the  power  entering  the  medium.  These  powers  can  be  calculated  quite  simply  for  a 
cylindrical  microwave  beam  of  radius  b  propagating  through  a  homogeneous  plasma  of 
length  £  as  depicted  in  Fig.  (A-l).  The  power  flowing  through  surface  S]  is 


(A- 18) 
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and  the  power  flowing  through  surface  Sj  is 


Figure  A-1.  Cylindrical  microwave  beam  with  radius  b 
traveling  through  a  homogeneous  plasma. 

where  the  subscripts  1  and  2  refer  to  the  value  of  the  electromagnetic  fields  at  surfaces 
S|  and  S2  respectively.  The  brackets  <  >  signify  time-averaged  values.  For  the  cylindrical 
geometry  depicted  in  Fig.  (A  I)  the  elements  of  surface  area  dS|  and  dSi  are  defined 
as 


■F  |  -  d>  1  rilhflr 

The  electric  and  magnetic  fields  at  S]  (i.e  .  z  =  0)  are 

1  1 0  1)  1  I  1  uj.  1 1 1  li 

1  ’  \  l> 


(A-20) 


( A-2 1 ) 


and 
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BifO.t)  =  i  —  [ B  uos  (wt)  -  a  sin  (cut)] 

1  >  (O 


Therefore, 


PK  2tt  1> 

-sr.U  rm' 


-  O 


(A-22) 


The  electric  and  magnetic  fields  at  S2  (i.e.,  z  =  2)  are 


K2<l\0  =  i  sFoe-af!  cos  (cut  - /3a 


(A-23) 


and 


p  ^at 

B2(P.t>  =  1,  ~  [,3  cos  (cut  -  fiV)  +  a  sin  (cut  —  /3I01 


Therefore. 


—  2af  „  , 

P'V'  277  h 

P  .  — — - -  /  /  rdSdr  =  ffb‘  .  . 

n"t  2cofi(>  o  a  \  2 <u/c0 


J  V  2cu«  / 


(A-24) 


The  attenuation  in  units  of  tlb  is 


Mdb)  =  -]0  log'°\T^)  =  20<,n°SlO  e  =  8.686  af  (A-25) 

It  should  be  emphasized  that  the  derivation  of  Eq.  (A-25)  is  based  upon  the 
propagation  of  a  plane  electromagnetic  wave  through  a  homogeneous  plasma  for  which 
Ohm’s  law  (J  =  oE)  is  applicable.  (Also,  boundary  effects  have  been  ignored  in  this 
derivation).  The  important  condition  here,  however,  is  homogeneousness  of  the  plasma 
(i.e..  Vcr  =  0).  This  aspect  is  often  overlooked  (see  for  example,  Ref.  16).  Thus  the  inversion 
of  the  total  attenuation  A(db)  into  an  absorption  coefficient  using  this  expression  is 
subordinate  to  these  conditions.  The  error  in  neglecting  these  conditions  has  yet  to  be 
evaluated. 
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NOMENCLATURE 

Anm  Transition  probability  for  transitions  from  the  nth  level  energy  state  to  the  mth 
level  energy  state 

A(x)  Microwave  attenuation,  lOlogio  (Poui/Pin) 

c  Speed  of  light  in  a  vacuum,  2.9979  x  108  m/sec 

Ej  Ionization  energy  level 

E„  Energy  of  the  nth  level 

e  Charge  on  electron  1.6021  x  10~19  coulombs 

G(y>  2  jTX1  f(r)  dx 

gn  Degeneracy  of  the  nth  energy  state 

h  Plank's  constant,  6.6256  x  10  27  erg-sec 

Inm  Intensity  of  the  radiation  from  nth  energy  level 

j 

mc  Mass  of  electron,  9.1091  x  10_31kg 

rig  Electron  number  density 

ng  Heavy  particle  number  density 

Pin  Microwave  input  power  to  plasma 

Pout  Microwave  power  at  exit  of  plasma 
Qen  Electron-neutral  collision  cross  section 

R  Radius  of  the  plasma 

r  Radial  coordinate 

Tc  Electron  temperature 

x,y  Cartesian  coordinates 
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Za +c  Electron  partition  function  for  the  ionized  atom 

Z^lec  Electron  partition  function  for  the  neutral  atom 

a(r)  Microwave  attenuation  coefficient 

enm  Emission  coefficient  for  spectral  emission  from  upper  energy  state  n  to  lower 
energy  state  m 

e0  Free-space  permittivity,  1 09/36jr  f/m 

A  1.23  x  107  Tc3/2  nc-1/2 

X  Wavelength 

Ho  Free-space  permeability,  4tt  x  10°  hy/m 

v  Total  collision  frequency,  vc{  +  ven 

vt j  Electron-ion  collision  frequency 

ven  Electron-neutral  collision  frequency 

a  Electrical  conductivity 

gj  Angular  frequency  of  microwave  signal 

cop  Angular  plasma  frequency,  (nee2/mee0 )% 
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